Order parameter oscillations in Fe/Ag/Bi 2 Sr 2 CaCu 2 08+5 tunnel junctions 
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We have performed temperature dependent tunneling conductance spectroscopy on 
Fe/Ag/Bi2Sr2CaCu208+s (BSCCO) planar junctions. The multilayered Fe counterelectrode 
was designed to probe the proximity region of the ab-plane of BSCCO. The spectra manifested 
a coherent oscillatory behavior with magnitude and sign dependent on the energy, decaying with 
increasing distance from the junction barrier, in conjunction with the theoretical predictions 
involving d-wave superconductors coupled with ferromagnets. The conductance oscillates in 
antiphase at E = and E = ±A. Spectral features characteristic to a broken time-reversal 
pairing symmetry are detected and they do not depend on the geometrical characteristics of the 
ferromagnetic film. 

PACS numbers: 74.45.+C, 74.50.+r, 75.70.-i 



In recent years numerous theoretical and experi- 
mental studies have been devoted to the spin polar- 
ized transport in ferromagnet(FM)-superconductor(SC) 
structures. This interest is equally motivated by the 
rapidly developing applications, e.g. in quantum two- 
level systems (qubits) based on phase shifts 0], as well 
as by the very interesting physics involving unconven- 
tional superconducting states. The exchange energy h 
in the ferromagnet leads to a spatially inhomogeneous 
state with the order parameter presenting phase changes 
by 7r in the proximity region at the interface with the 
SC. Such a modulated state is called Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) from the names of the authors who 
first predicted it |2(. In the FM, the Cooper paired elec- 
trons have different energies and Fermi momenta, so they 
dephase on a scale of several ferromagnetic coherence 
lengths. Consequently, this region presents a diversity 
of intensely explored phenomena like spatially dumped 
oscillations of the density of states (DOS) and gapless su- 
perconductivity oscillations of the superconduct- 
ing critical temperature T c and critical current I c with 
the thickness of the FM electrode in layered systems 
0iQ, or asymmetric zero- bias conductance peak (ZBCP) 
splitting 
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An extensive attention was dedicated to the LOFF 
state arising in junctions between conventional SC's and 
FM films of variable thickness. As example, the measure- 
ments on Al/A^Oa/PdNi/Nb junctions reported by the 
Orsay group clearly show the difference in the conduc- 
tance spectra at two different locations in the proximity 
region, indicating the order parameter change from pos- 
itive (0-state) to negative values (estate). However, the 
data originated from experiments involving high-Tc su- 
perconductors (HTSC) is rather scarce, even if it could be 
conclusive about the interaction between ferromagnetism 
and a more complex form of superconductivity. In partic- 
ular, the superconductivity mechanism in the cuprates is 
far from being fully understood. According to a series of 
phase sensitive experiments [lo| the symmetry of the or- 



der parameter is d-wave, pure or mixed with a secondary 
s or d component [Tlj . It also seems that the antifer- 
romagnetic spi n plays an important role in the pairing 
mechanism 12]. Therefore, it is constructive to consider 
that the LOFF state in the cuprates, while intrinsically 
interesting, may also offer information about what makes 
the high-Tc superconductivity different from the conven- 
tional one. A multitude of theoretical models were devel- 
oped to calculate the specific transport properties, spin 
and charge relaxation in tunneling experiments involv- 
ing d-wave HTSC with FM counterelectrodes. The most 
popular ones are those based on the Blonder-Tinkham- 
Klapwijk (BTK) approach [H[li|. However, these mod- 
els consider just the transmission and reflection processes 
at the interface with the FM film modeled as a semi- 
infinite domain. In this paper we will use a model || 
based on a ballistic quasiclassical formalism with the 
DOS averaged over all classical trajectories of multiple 
reflections in the FM domain of finite thickness d. 

Hoping to contribute to the experimental information 
about this topic, we report in this paper the transport 
properties in a multilayered-Fe/Ag/Bi2Sr2CaCu208_<5 
(BSCCO) ab-plane oriented planar tunneling junction. 
The BSCCO ab-plane is probed successively from 5 over- 
lapping Fe thin films with the same thickness, through 
an intermediate thin layer of Ag designed to minimize 
the barrier strength and avoid the appearance of a spin- 
glass phase at the junction interface |l5l [. The spectra 
are collected on the thicker and thicker counterelectrodes 
and studied as function of the temperature and thick- 
ness of the FM domain. The planar junctions were built 
on freshly cleaved monocrystal samples of slightly un- 
derdoped BSCCO grown by the self-flux method. Two 
silver paste leads were attached to the crystal and its 
lower side was molded in epoxy resin. Before the epoxy 
dried, the upper side was fixed between two quartz glass 
slabs with coplanar surfaces. The epoxy penetrated be- 
tween the sides of the slabs so that the crystal was 
caught between two smooth glass surfaces attached by 
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FIG. 1: Junction layout. The layered Fe/Ag conterelectrode 
is deposited perpendicular on the BSCCO ab-plane. 



thin epoxy walls. While the surface of the glass was 
protected, we performed a clean mechanical polishing 
along the ab-plane of the crystal, laterally embedded in 
epoxy. We monitored the process with a profilometer, 
until the exposed edge of the crystal was at the same 
level as the glass surfaces. The counterelectrode sand- 
wich was then evaporated using a metallic shadow mask 
shaped as in Fig. ^ forming thin film strips 0.1 mm in 
width. The deposition was controlled using a calibrated 
quartz thickness-monitor. Above a 30A silver film, five 
30A iron films were deposited at a rate of O.lA/s for uni- 
formity [lfij . The residual resistivity ratio of the Fe film 
RRR = R(300K)/R(A0K) « 8.3 indicates a good de- 
position quality. The tunneling barrier is formed natu- 
rally on the exposed edge of BSCCO. It had a resistance 
Rb(300K) 450S1, for the typical junction presented 
here, showing a high transparency. We evaluated the 
elastic mean free path I of conduction electrons in one 
of our Fe films using an expression derived from Pippard 
relations 0: v F l = (irk B /e) 2 (cr/ 7 ). With the iron film 
conductance a «4xl0 3 S/m 2 , electronic specific heat 7= 
5xl0~ 3 J/K 2 , and Fermi velocity vp = 1.98xl0 6 m/s, we 
obtained I wlOOA. Since the quasiclassical description of 
conductance assumes I < d 0| , we see that such a model 
is applicable to our first three Fe layers. We performed 
conventional 4-point measurements, with a constant cur- 
rent driving bias. 

The tunneling spectra at the interface with the BSCCO 
crystal were obtained successively for the six layers, start- 
ing with the Ag counterelectrode, then probing an in- 
creasing thickness in the FM region from to 150A, by 
30Aincrement for each electrode. Each layer was mea- 
sured at different temperatures, from 150 to 4.2 K. Fig. 2 
represents the temperature dependence for the first three 
layers, each curve being normalized by the spectrum at 
150K. All layers manifested a sudden change in the back- 
ground conductance at T c . This phenomenon can be 
attributed to the competition between a screening spin- 
accumulation close to the interface, responsible for a drop 
in the barrier conductance, and the onset of Andreev 
reflections, with an adverse effect 01 ■ I R an Andreev 
reflection, the electron incident on the barrier with en- 
ergy lower than the superconductive gap is reflected as a 




FIG. 2: Temperature dependence of the conductance spectra. 
The arrows indicate the background conductance variation at 
T c as the temperature is decreased, a) Normal layer b) 3uA 
iron layer, and c) 60A iron layer. 



hole, while a Cooper pair is transmitted into the SC. In 
the case of d-wave superconductors the Andreev reflec- 
tions are phase sensitive. This leads to the formation of 
bound states at the junction interface with amplitudes 
dependent on the angle of incidence. In a quasiclassical 
picture 5], the particles move along trajectories reflected 
between the boundaries of the counterelectrode, the den- 
sity of states per trajectory being correlated to the length 
of the trajectory or number of reflections which in turn 
depends on the energy and exchange potential. The tun- 
neling conductance spectrum is a measure of the energy 
resolved total density of states and is regarded as a sum 
of local bound states averaged over all trajectories dis- 
tributed by length. Therefore, the thickness of the coun- 
terelectrode influences the amplitude of the bound states 
and thus the spectral features, like the ZBCP. For our Ag 
terminal, since the amplitude of the Andreev states is low 
(Fig. [21 and |2J), and it is in direct contact with the next 
iron layer, the background conductance at T c decreases 
with decreasing temperature, due to the predominance 
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FIG. 3: Tunneling spectra nor- 
malized by the background conduc- 
tance Go, represented for increas- 
ing counterelectrode thickness d, at 
three temperatures T. a) T = 5 K. 
The arrow on the Ag/BSCCO spec- 
trum points to the minigap edge. 
The asymmetric spectrum at d = 
60A permits an evaluation of the Fe 
film polarization from the peak ra- 
tio, P Fe « 43%. b) T = 40 K spec- 
tra (continuous lines) compared to 
the theoretical fits (dash lines) with 
h = 170 meV (20 meV for the first 
Fe layer), A = 44 meV, and Z = 
0.51. The fitted d values (Fig. 
match close enough the actual val- 
ues, c) T = 60 K. The ZBCP for d = 
150A is not split at this higher tem- 
perature approaching T c , due to the 
suppression of the BTRS and the big- 
ger effect from the spin-splitting field 
on the superconductivity since A de- 
creases while the h remains constant. 



of a spin accumulation at the junction barrier. The tem- 
perature dependence of the Ag/BSCCO tunneling spec- 
trum shows no ZBCP (Fig. This may indicate a 
predominant injection normal on the [100] plane (ll^Hif 
combined with a rarefied distribution of long trajectories 
due to the very low thickness of the film, leading to the 
suppression of low energy states. The energy resolved 
DOS shows signs of a minigap (Fig. |3J) at zero energy 
bias, which points out to some degree of disorder in the 
counterelectrode favorable to short trajectories or higher 
energy subgap bound states. A pseudogap opens at T* rj 
160 K and the superconducting transition takes place at 
T c k 85 K, clearly visible on the T-V plane contour. We 
measured a peak to peak SC energy gap A ss 44 meV at 
40 K. 

Next we probed the FM region of the junction. With 
each Fe layer, 30A are added to the total distance ac- 
cessible to the pairs injected into the counterelectrode. 
The pairs survive along a distance of several coherence 
lengths £ w 38A as the SC order parameter decays oscil- 
latory around zero. We performed theoretical fits for the 
DOS spectra at the intermediate temperature 40 K. The 
model we employed Q calculates the DOS taking into ac- 
count the thickness d of the FM film when the interface 
with the SC is highly transparent and d is smaller than 
the electron mean free path I. With I calculated above, 
the model is applicable mostly for the first few layers, 
which are anyway the most suggestive. Since we expect 
our boundaries to be rough, leading to an admixture of 
trajectories and a smearing of the spectral features, the 
DOS was averaged over a Gaussian distribution of thick- 
nesses around each d. The best fits imposed a value h 



= 20 meV for the first iron layer and h = 170 meV for 
the second one (consistent with the accepted value [l6|l- 
The fits for all other Fe layers were obtained by keeping 
constant the exchange energy h = 170 meV, the gap A = 
44 meV, and the barrier transparency Z = 0.51. The only 
adjusted variable was the thickness scaled by the coher- 
ence length, d/n^p. The procedure resulted into five val- 
ues for the layer thickness remarkably close to the actual 
values (see Fig. QJd) , confirming the validity of the model 
and also the regularity of the thin films. We also take 
advantage from the BTK type models to underline the 
effect of the d-wave (possibly mixed with s-wave) pairing 
symmetry of BSCCO. 

The temperature dependent spectra measured on the 
first FM layer (d — 30A) (Fig. also has a pseudogap. 
The T c is slightly smaller and the background conduc- 
tance after the transition increases, indicating the grow- 
ing influence of the Andreev reflections. The most re- 
markable features (Fig. |3J) , are the dips replacing the ± A 
coherence peaks. Due to the direct contact with the Ag 
layer, the spectrum rather shapes up the effect of a distri- 
bution of spin-polarized quasiparticles in the Ag film, so 
that the theoretical fit on this curve imposed an exchange 
energy h=20 meV, lower than that for the other Fe layers. 
For 0< h < A, the ±A peaks move to subgap energies, 
being replaced by minima in the DOS |5(. The resulted 
subgap peaks are expected to merge into a ZBCP with 
increasing d, and to exhibit local oscillatory variations. 
Indeed, the next FM layers show the oscillations, but the 
peaks do not merge with increasing d. This may be due 
to properties of the junction unaffected by the thickness 
of the FM environment, like the presence of a mixed su- 
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FIG. 4: The reduced DOS develops coherent oscillations de- 
pendent on d. a) The experimental results approach better 
the predicted behavior for increasing temperature, due to the 
decreased s-wave component, b) There is a remarkable match 
between the theoretical thicknesses and the actual d values. 



perconductive phase with broken time-reversal symmetry 
(BTRS) at the junction interface (e.g., (d+is)-wave) 111. 
This splits the ZBCP, forming a gaplike feature with the 
size given by the magnitude of the sub-component s, and 
influenced by temperature. In this case, the s-wave gap 
remains unchanged with increasing d, while the subgap 
DOS oscillates up and down as the decaying order pa- 
rameter alternatively takes negative and positive values 
respectively. 

As seen in Fig. the next Fe layer (d = 60A) mea- 
sures a natter pseudogap at high temperatures. The con- 
ductance enhancement at the superconductive transition 
is larger, consistent with the high amplitude of the An- 
dreev peaks. The electronlike and holelike quasiparti- 
cle branches of the spectrum are asymmetric at lower 
temperatures. The split ZBCP is high and it becomes 
broader and less asymmetric with increasing tempera- 
ture (Fig. |2^,c). This asymmetry may be attributed 
to the presence of a spin dependent contribution to the 
ZBCP split at lower temperatures. If a spin imbalance 
is present in the barrier, the ZBCP corresponding to up 
(down) spin currents are shifted at higher (lower) ener- 
gies. The polarization in the FM layer decreases the spin- 
up peak and increases the spin-down peak. The ratio of 
the peaks can be used to evaluate the spin polarization P 
in iron. In our case, Pp e ~ 43%, which is consistent with 
the values reported in literature 0|. The broader ZBCP 
at temperatures closer to T c may be the effect of a ther- 
mal smearing and the suppression of superconductivity 
due to the diminishing A. 

The next three layers (d = 90-150A) complete an os- 
cillatory cycle. The curves are too noisy to discern the 
more delicate features but the general behavior is suffi- 



ciently clear: the superconductivity is highly suppresed 
and the LOFF states develop d-dependent coherent os- 
cillations, with smaller and smaller amplitude. The sign 
of these oscillations depend on the energy. In Fig.^J the 
reduced DOS, G(E)/Gq — 1, oscillates in antiphase at 
energy E = with respect to the E = ±A values. At 
lower temperatures the zero energy conductance is lower 
than the theoretically predicted curve (Fig. 0^,), due to 
the significant presence of the BTRS induced ZBCP split. 
The BTRS states disappear at a temperature lower than 
BSCCO T c , so that the 150 A spectrum presents a 60K 
ZBCP without the split. Consequently, the points taken 
at temperatures nearer to T c approach better the theo- 
retical expectation. 

In summary, we constructed planar Fe/Ag/BSCCO 
junctions, with the Fe region imparted in layers of equal 
thickness. Temperature dependent tunneling spectra 
were collected on each layer, probing a reproducible spa- 
tial dependent instance of the LOFF state. The spectral 
features are affected by the pair-braking exchange field 
in the FM region, the thickness of this area and the un- 
conventional pairing symmetry in the BSCCO. As the 
Fe thickness is increased, the conductance spectra show 
decaying oscillations with energy dependent amplitude 
and sign. The ZBCP is split regardless the thickness, 
but dependent on the temperature, due to the presence 
of a subcomponent added to the dominant d-wave com- 
ponent at the junction interface to form a BTRS state. 
The spectrum is asymmetric at low temperature due to 
the spin-polarization. 

This work is supported by NSF Grant No. 
DMR9972071. 
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